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Abstract 

r) ' A method of controlling Shil'nikov's type chaos using windows that appear in the 

^ . 1 dimensional bifurcation diagram when perturbations are applied, and using exis- 

tence of stable homoclinic orbits near the unstable one is presented and applied to 
the electronic Chua's circuit. A demonstration of the chaos control in the electronic 
circuit experiments and their simulations and bifurcation analyses are given. 
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1 Introduction 



Controlling cliaos is important for application in the comunication systeni[l,2] 
and controlling dynamics in technology like convective flows, laser excitation 
and electronic currents. There are two well-known methods for controlling 
chaos 1) OGY method[3] and 2) feedback (Pyragas) method[4]. In the case of 
electronic current, chaotic oscillation of Shil'nikov type[5] occurs in the Chua's 
circuit and the chaotic patterns depending on the parameters of the dynamical 
system were extensively investigated [6]. 

Chua's circuit consists of a autonomous circuit which contains three-segment 
piecewise-linear resistor, two capacitors, one inductor and a variable resistor. 
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The equation of the circuit is described by 



^^di^ci = G{vc2 - vci) - g{vci,mo, mi] 

C2J1VC2 = G{vci - VC2) + IL 



LiiL = -VC2 



where vci and Vc2 are the vohages of the two capacitors (in V), il is the 
current that flows in the inductor (in A), Ci and C2 are capacitance (in F), L 
is inductance (in H) and G is the conductance of the variable resistor (in VL~^). 
The three-segment piecewise-linear resistor which constitutes the non-linear 
element is characterized by 

g{vci,mQ,mi) = (mi - mo){\vci + Bp\ - \vci - Bp\)/2 + rrioVci, 

where Bp is chosen to be IV, mo is the slope (mA/V) outside |fci| > Bp and 
mi is the slope inside —Bp < vci < Bp. The function ^(fci,mo,mi) can be 
regarded as an active resistor. If it is locally passive the circuit is tame but 
when it is locally active it keeps supplying power to the external circuit [13]. 
The chaotic behavior is expected to be due to the power dissipated in the 
passive element. A difference from the van-der Pol circuit is that the passive 
area and the active area are not separated by the strange attractor, but they 
are entangled. 

We control chaotic circuits by triggering simple electronic pulse as [8]. Main 
difference from this work is that we choose a proper height of the pulse, so 
that the periodic trajectory is kept and we do not need to reset the current. 

By coupling two Chua's circuit, one can control the chaotic motion by phase 
synchronization [9]. We adopt the phase synchronization technique on a single 
Chua's oscillator and control the chaotic motion. We adopt also the feedback 
control technique to the Chua's oscillator. 

The Chua's electronic circuit and connection of the comparator and the pulse 
generator are shown in Fig.l. Changing the scales of variables as Ci x 10^, L x 
20, G X 10^, mj X 10''^ and t x 2 x 10"^, we obtain the differential equation 

—X = f ■ X where x = (x(r), |/(r), z{t)) as 



—x = a{y-x- g{x, mo, mi)) 



£y 



dr' 



X — y + z 



(2) 
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Fig. 1. Chua's electronic circuit with impulse generator. 

C2 



and 



g{x, mo, mi) = (mi — mo){\x + 1| — |a; — l|)/2 + rriox. 
We choose Ci = |, C2 = 1, L = f, G = 0.7, mo = -0.5, mi = -0.8, B„ = l. 



;, -t^p 



The differential equation is invariant under symmetry {x, y, z) -^ {—x, —y, —z) 

and the state space can be divided into three domains by the planes f/_i = 

{z = —1} and f/+ = {z = 1}. In each domain there are equilibrium points 

mo — mi ^ 
at P+ = {~k,0,k) O=(0,0,0) P^{k,0,-k) where k = —-. In the 

TTlg + i 

three regions, the linearized f has a real eigenvalue 7 and a pair of complex 
eigenvales a ±iu. 

Shil'nikov considered a continuous vector field X on R^ with equilibrium 
points p such that \a\ < \'j\ and u ^ and that there is a homoclinic or- 
bit from p. Then he proved that there is a perturbation Y oi X such that Y 
has invariant sets containing transversal homoclinic orbits[10,5]. It means if 
the Chua's circuit satisfies the condition \a\ < \'~f\ the trajectory can be per- 
turbed into an orbit with countable sets of horseshoe, or Chaos is embedded 
in the system. 

Newhouse showed, in addition to bifurcations mechanism of Shil'nikov, subtle 
and complicated dynamical behaviour associated with homoclinic tangency 
exists in the diffeomorphism of continuous vector field X on R^. He defined a 
zero-dimensional hyperbolic invariant set A in a plane, stable manifold W''^{A), 
unstable manifolds iy(A), thickness of W^*(A) as r*(A) and that of Vr"(A) 



as r"(A). He showed when 1) r*(A2) ■ r"(Ai) > 1, 2) W'{K2) n l^"(Ai) and 
W^"(A2) n W^*(Ai) both have points of transversal intersection which are not 
in Ai U A2, 3)W^''(A2) and W^"(Ai) have a point of tangency, then the diffeo- 
morphism has infinite number of sinks[7,10]. 

The Shil'nikov's chaos occurs in the R^ but the theorem on 2-dimensional(2-d) 
Cantor set would apply also to 3-d Cantor set since the homoclinicity of the 
orbit is the same. 

A sufficient condition for finding a map f in the neighborhood of p such that for 
m > M, f™" has an invariant set A^ topologically equivalent to the horseshoe 
is that the eigenvalues of the Jacobi matrix at the fixed point p and A satisfy 
p < 1 < A < p-^[10]. When r*(A2) ■ r"(Ai) > 1 is satisfied and 1^"(A) has 
a point of tangential intersection with W^"(A), A is called wild hyperbolic set 
and the region around this set is called Newhouse region. 

In the windows region of the 1-d bifurcation diagram as a function of the 
strength of the perturbation, there are countable periodic orbits belonging to 
W^*(A2) and to W^"(Ai). Since it is plausible that there are points of transver- 
sal intersection outside Ai U A2 in the windows region, perturbation in this 
region would kick the unstable orbit to one of the sinks. We define the strength 
of the pulse or the amplitude of synchronizing oscillation or the amplitude of 
the feedback by using the information of the position of windows in the 1-d 
mapping of the bifurcation diagram. Various kinds of shift from unstable man- 
ifold to stable manifold occurs due to the dense presence of stable homoclinic 
orbits near each unstable orbits in the Newhouse region[10]. 

This paper is organized as follows. In sect. 2 we explain the experimental setup 
of the Chua's circuit control system and chaotic patterns. The method of 
controlling the double scroll and the experimental results are shown in sect. 3. 
In sect. 4 we show the method of bifurcation analysis. Impulse control of the 
double scroll is analyzed in sect. 5. Drive control and feedback control are 
analyzed in sect. 6 and in sect. 7, respectively. Analysis of the spiral control is 
given in sect. 8. Discussion and outlook are presented in sect. 9. 



2 The Chua's circuit 



The dynamical system of the Chua's circuit is defined as eq.(2). When the 
variable resistance is about R = 1616^7, on the oscilloscope displaying vci — 
Vc2, a circle appears and as the resistence is reduced to -R = 1534^2 suddenly 
a chaotic spiral as shown in Fig. 2 appears. When R is reduced to R = 1528fi 
a 3 cycle spiral appears and when R is 1457^2 a double scroll as shown in 
Fig. 3 appears. When R is further reduced to 1402^2, a periodic double scroll 



Table 1 

Eigenvalues of the Jacob! matrix at the fixed point. 
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with 3 cycles around each fixed point appears and for smaller R the double 
scroll disappears. By measuring Lyapunov exponents we check that the spiral 
and the double scroll are chaotic, and the chaotic behavior is implied by the 
presence of the period-3 cycle according to the Li-Yorke's theorem[10]. 



df 



(x,y,z) for \x\ < 1 are 



Eigenvalues of the Jacobi matrix 7— where x 

ax 
defined as o"o =t ?cJo,7o and for x < —1,1 < x are defined as o"i ± iui,-yi. 

Numerical values as a function of the resistence R used in the experiment and 

the conductance G used in the simulation are shown in Table 1. The value 

G used in the simulation does not necessarily coincide with the experimental 

value 1000/i?. The Shil'nikov's condition for the homoclinic chaos is \cr\ < \'~f\ 

and cj 7^ and the condition is satisfied in the range of 0.55 < G < 0.643. 



The chaotic orbits can be guided to periodic orbits by triggering electric pulse 
as shown in Fig. 4 for the spiral and Fig. 5 for the double scroll. In these photos, 
the abscissa is vci and the ordinate is vc'2- The Figs. 6 and 7 are the time series 
of the vci and the Vc2, respectively. 



Since the electric field pulse applied to the system is not a delta function but 
contains a width, the shift 6y is not exactly equal to the height of the pulse. In 
the case of double scroll, the number of cycles in z > region and z < region 
depend on the voltage Vc2 of the sink that the orbit from unstable manifold 
is trapped, and so the exact quantitative simulation is impossible. But we can 
study qualitative mechanism of the chaos control by simulation. 





Fig. 2. The Chua's spiral. 



Fig. 3. The Chua's double scroll. 





Fig. 4. The Chua's spiral impulse-1 Fig. 5. The Chua's double scroll im- 

type control. pulse-1 type control. 





Fig. 6. The voltage of the capacitor-1 
in double scroll impulse-l type control. 

3 Controlling the double scroll 



Fig. 7. The voltage of the capacitor-2 
in double scroll impulse-1 type control. 



The Poincare map on the plane of x = 1 of the double scroll orbit obtained 
by using parameters after eq.(2) is shown in Fig.8. The points in the inner 
circle (region a : —1 < z < —0.4) move to the outer branch (region h : 
—0.3 < y < -0.2,-1.3 < z < —0.7) and go to the top of the outer circle 
(region c : 0.1 < y < 0.2) and go to region d or to its left neighbor (region 



e : —0.1 < y < 0.1) and go to complicated routes. When the trajectory goes to 
the latter region, the trajectory becomes complicated. Thus in order to realize 
periodic orbit, it is necessary to avoid the points trapped into region e. 

The double scroll pattern of the Chua's circuit can be controlled by triggering 
electric pulses to the system. We show two types of triggering, i.e. one pulse 
as the trajectory crosses x = 1 (impulse- 1) and pulses as the trajectory passes 
X = 1 and X = —1 (impulse-2). 

When one triggers the impulse 5y when the orbit crosses x = 1 from x larger 
than 1, the curvature radius of the trajectory becomes larger. The points that 
the route takes for 5y = —0.041 are points A, and the points for 5y = —0.06 
are the spiral that passes through B, the points for Sy = —0.08 are the spiral 
that pases through C, and the points for 6y = —0.12 are the spiral that passes 
through D. Isolated points E, near {y,z) = (—0.6,-1.8), (-0.3,-0.7) and (0,- 
0.7) are for 6 = —0.14. The points that are trapped in —0.1 < y < 0.1 region 
become shifted to 0.1 < y < 0.2 region, as shown in Fig.9. 



/ 



I i 



/ 1 



I 



-o.i-c.a-a,2-D.i D ■:.! 0.2 
y 



:■: 


^ 


< 


r\ 


'" 


.J 


• !••• 


f « 


• 


/" 




DCB 

> 




4^ 





-O.G -0.4 -0.2 0.2 

y 



Fig. 8. The Poincare map of the double Fig. 9. The Poincare map of the double 
scroll. scroll after triggering the impulse. 



Using the above information, we trigger an electric impulse when the orbit 
passes the plane x = 1, and study 1-d bifurcation diagram as a function of the 
height of the impulse. 
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Fig. 10. The 1-d bifurcation diagram as 
a function of the impulse strength 5y of 
impulse-1 type. 



Fig. 11. The 1-d bifurcation diagram as 
a function of the impulse strength 6y of 
impulse-2 type. 



We add a shift 6y as the trajectory passes the x = 1 plane, and call it the 
impulse- 1 type. In a plot of x coordinate of the y = section as a function of 
the Sy, we find windows for 6y = —0.41, -0.09 and -0.14 as shown in Fig. 10. 
In the impulse-2 type, we add a shift 6y as the orbit passes the plane x = 1 
and —6y as the orbit passes the plane x = —1. The windows of the impulse-2 
type are shown in Fig. 11. 

The appearance of windows suggests a presence of heterocrinic orbits[17], or 
change of unstable orbits to stable orbits in a certain parameter region. 



4 The bifurcation analysis 



Chaotic system changes its Poincare map when certain parameter of the dy- 
namical equation is changed, and this phenomenon is known as bifurcation. 
In the present Chua's circuit, the parameters a and P are fixed, but the 
strength of the additional electric pulse is a parameter that changes a com- 
plicated orbit to a torus. In order to analyze this qualitative changes, we 
consider the Poincare mapping Tx i.e. a diffeomorphism from the state space 
u = (mi, U2, Ms) into R^ 



<^ 



n,(P2,(P3) 



Tx{u) 



(3) 



As in [12], we consider the first derivative of the Poincare mapping 



|.r.(u) = BT.(u) = |^(..u,A) 



(4) 



which is a tensor 



d(f) 
du 

that satisfies the diff^erential eq. 
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and 



dg_ 
dx 



(5) 



(6) 



rrii 
mo 



X < -1 
-1 <a; < 1 
X > 1 



(7) 



The differential equation eq. (6) is solved by the 4th order Runge-Kutta method 

d(t) 
with initial condition — — = / at initial r to r + Tr,er, where r is chosen at the 

du 
point when the pulse in y direction is triggered and Tp^r is the period between 

the pulses are triggered. We use data u(r) of the periodic trajectory for fixing 

the time that the pulse is triggered and thus the y coordinate is shifted. We 

measure eigenvalues of — — at each r and store. 

du 

In most cases, after a few hundred steps, one of the three eigenvalues of the 

first derivative of the Poincare mapping damps to very close to and the 

(9f 
diffeomorphism defined by 7— becomes essentially 2-dimensional. 

ax 



5 Impulse control of the double scroll 



We simulate the Shil'nikov's chaos by choosing the initial condition u(l) = 
(0.01,0.01.0.01) and take dr = 0.01 and make 50000 steps. Using the param- 
eter a = 9, P = 14.28 we obtained the double scroll Fig. 12. 

We consider 2 types of controlling Shil'nikov's double scroll by impulse. 

• Impulse-la type: When x passes 1 from x„ > Xn+i shift y to y — 0.041. 

• Impulse-lb type: When x passes 1 from x„ > Xn+i shift y to y — 0.097. 

• Impulse-lc type: When x passes 1 from x„ > Xn+i shift y to y — 0.147. 

• Impulse-2a type: 

a) When x passes 1 from Xn > Xn+i and y < —0.37 shift y to y — 0.05 



b) When x passes —1 from x„ < x„+i and y > 0.37 shift y to y + 0.05 
Impulse-2b type: 

a) When x passes 1 from x„ > x„+i and y < —0.37 shift y to y — 0.18 

b) When x passes —1 from Xn < a;„+i and y > 0.37 shift y to y + 0.18 



5.1 The impulse-la type control 



The simulation results of the controlled trajectories (the last 5000 steps of the 
50000 steps) by the impulse- la type is shown in Fig. 13. 





Fig. 12. The trajectory of the double Fig. 13. The trajectory of the double 
scroll without impulse. scroll controlled by impulse-la type. 



The X coordinate of the trajectory obtained by applying the shift 6y 
are shown in Figs. 14 . 
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The first derivative of the Poincare mapping is calculated from r = 1838 
when the Sy is applied to r = 3705 when the subsequent 6y is applied. Its 
eigenvalues consist of a pair of complex and one real. The behaviors of the 
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real eigenvalue and two complex eigenvalues before and after the second pulse 
are shown in Fig. 15 and in 16. The imaginary part of the complex eigenvalues 
crosses Im/i = ±1 near r when the impulse is applied. The real part of the 
complex eigenvalue is positive but the slope of increasing tendency is reduced 
when the impulse is applied. 
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Fig. 16. Magnification of the small- 
est eigenvalue of the first derivative of 
the Poincare mapping of the impulse-la 
type. 



Fig. 17. Tangency of the controlled or- 
bit (impulse-la) and the homoclinic or- 
bit. 



The trajectory of the impulse-la type suggests that the orbit in the unstable 
manifold W"^ is perturbed in the region A2 and shifted to the stable manifold 
W^ and near the fixed point P+ = (—1.5,0,1.5) crossing to the unstable 
manifold W"^ occurs in the region Ai. In Fig. 17, we show the homoclinic orbit 
started from the vicinity of the fixed point P| = (—1.5 + 0.001, 0, 1.5) and the 
controlled trajectory of the impulse-la type. The figure suggests that after the 
impulse, the trajectory is absorbed into the homoclinic orbit and the tangency 
of the unstable manifold PF"(A_) and the stable manifold PF*(A+) occurs, 
where A± implies the region near the fixed point P±, respectively. Products of 
the three eigenvalues after r = 2035 are almost constant of 0.0069(1), which 
implies that the orbit is on W^{A^) It is interesting that the absolute value of 
the imaginary part of the eigenvalue becomes smaller than 1 as the trajectory 
shifts to the manifold PF'^(A+). 



5.2 The impulse- lb type control 



The second window in the 1-d bifurcation diagram exists near 6y = —0.097. 
The X coordinate of the trajectory for this perturbation is shown in Fig. 19. 
There are 3 cycle oscillation in y > and 3 cycle oscillation in y < 0. 

The behavior of the eigenvalues of the first derivative of the Poincare mapping 
in the case of 6y = —0.097 is different from that of —0.041. We measure it 
from r = 319 when the first 6y is apphed, to r = 1599 when the subsequent 
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Fig. 18. The trajectory of the Chua's cnxuit impulse-lb type. 
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Fig. 19. The x coordinate of the trajec- 
tory of impulse-lb type. 



Fig. 20. The real part of the eigenvalue 
of the first derivative of the Poincare 
mapping of the impulse-lb type. 



6y is applied. In Fig. 20, real eigenvalues of the first derivative of the Poincare 
mapping near the time when the second pulse is applied are shown. A crossing 
of the eigenvalue with the line of /x = 1 occurs as the pulse is triggered. One of 
the three eigenvalues of the first derivative of the Poincare mapping is zero and 
the product of the rest two eigenvalues before r = 1599 is constant (-0.170(2)), 
which means that the orbit has been on the iy(A_|_). 



5.3 The impulsel-c type control 



Corresponding to the relatively large window in the 1-d bifurcation diagram 
near Sy = —0.14, we trigger 6y = —0.14 in the impulse-lc type. There are 2 
cycle oscillation in y > and 2 cycle oscillation in y < 0. The x coordinate of 
the trajectory for impulse-lc are shown in Fig. 22. 

The first derivative of the Poincare mapping is calculated from r = 685 where 
the Sy is applied to r = 1609 where the subsequent 6y is applied. In Figs. 23, 
the eigenvalues are shown. Before the impulse, there is a pair of complex 
eigenvalues and after the impulse, a pair of real eigenvalues appear. Product 
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Fig. 21. The trajectory of the Chua's circuit impulse-lc type. 
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Fig. 22. The x coordinate of the trajec- 
tory of impulse-lc type. 
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Fig. 23. The real part ( branching at 
r = 1554) and the imaginary part (van- 
ishing at r = 1554) of the eigenvalue of 
the first derivative of the Poincare map- 
ping of impulse-lc type. 



of the three eigenvalues is constant (0.58(1)) before r = 1554 and the orbit has 
been on W^*(A_|_). After changing from complex eigenvalue to real eigenvalue, 
one eigenvalue becomes larger than 1 and become unstable. 



In order to check whether the transition at r = 1554 is smooth, we calculate 
the tetrahedron Af produced by mapping of 4 points (x — 0.1,?/, 2;), (x, y — 

di 

0.01, 2;), (x + 0.1, y, z), (x, y + 0.01, z) by the Jacobi matrix — at r = 1546, 

1553 to 1557, 1560 and 1564. ( Fig.24 ) The area projected on the x — y plane 
is nearly parallelogram . One of the eigenvalues of the first derivative of the 
Poincare mapping is consistent with and the rest are a pair of complex value 
in 1546 < r < 1553 and at r = 1554 all the eigenvalues become real. In the z 
direction, a shift of lS.z = —2 occurs between r = 1555 and 1556, but on the 
projected x — y plane the movement is smooth, which suggests that the orbit 
in W'^{P^) is smoothly connected in W^{PS). 



13 




Fig. 24. Mapping of the tetrahedron spanned by 4 points 
(x — 0.1, y, z), (x,y — 0.01, z), {x + 0.1, y, z), (x,y + 0.01, z) by the first derivative of 
the Poincare mapping for r = 1546,1553 to 1557, 1560 and 1564 from right to left. 
A shift of Az = —2 occurs between r = 1555 and 1556. 

5.4 The impulse-2a type control 



The iinpulse-2 type control is done by applying impulse at a; = —1 as well as 
X = 1. In this case, eigenvalues of the first derivative of the Poincare mapping 
are real near the point where the 2nd impulse is applied to the system. The 
real largest eigenvalue crosses /z = 1 from above near the time when the 2nd 
pulse is applied. The slope of the negative eigenvalue becomes flat at the same 
time. The 1-d bifurcation diagram shown in Fig. 11, we find 4 windows, but 
we present results of 2 windows. 

In the impulse 2a type control, we choose 5y = 0.05, whose controlled trajec- 
tory (the last 5000 steps of 50000 steps) is shown in Fig. 25. 

The first derivative of the Poincare mapping is calculated from r = 771 when 
the 6y is applied to r = 1927 when the subsequent 6y is applied. In this 
region one of the three eigenvalues is consistent with and the product of the 
rest two eigenvalues before r = 1927 is smaller than -1.5 but after r = 1927 
it rapidly increases and becomes about 3 at r = 1937. The tangency of the 
unstable manifold W^{A_) and the stable manifold VI^*(A_|_) can be visuahzed 
by comparing the controlled trajectory in the y — z plane and the homoclinic 
orbits started from the vicinity of the fixed point P^ = (—1.5 -|- 0.001, 0, 1.5) 
and P!_ = (1.5 — 0.001, 0, —1.5), which are shown in Fig.28 
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Fig. 25. The trajectory of the Chua's circuit impulse-2a type. 
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Fig. 27. Eigenvalues of the first deriva- 
Fig. 26. The x coordinate of the trajec- five of the Poincare mapping of the im- 
tory of impulse-2a type. pulse-2a type. 




Fig. 28. Tangency of the controhed orbit (impulse-2a) and the two homocHnic orbits. 
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5.5 The impulse-2h type control 



In the inipulse-2b type we choose 5y = 0.18. The trajectory is shown in Fig. 29 
and the projection on the a;-axis is shown in Fig.30. 
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Fig. 29. The trajectory of Chua's circuit impulse-2b type. 
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Fig. 30. The x coordinate of the trajec- 
tory of impulse-2b type. 



Fig. 31. Eigenvalues of the first deriva- 
tive of the Poincare mapping of the im- 
pulse-2b type. 



The first derivative of the Poincare mapping is calculated from r = 191 when 
the 6y is applied to r = 748 when the subsequent 6y is applied. Its eigenvalues 
are shown in Fig. 31. Product of two real eigenvalues are negative and smaller 
than -1.3 before r = 748 but it rapidly approaches to -1 after r = 748. 
Products of the three eigenvalues after r = 748 is constant (0.38(1) x 10~^), 
which means that the orbit is absorbed in the stable manifold PF^(P+). 



6 Drive control of the double scroll 



As in the forced oscillator, the chaotic system can be controlled by applying 
sinusoidal perturbation F sin lot. In this case in the next feedback control, we 
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choose G = 0.65 so that the double scroll is produced. 



£x = a{y - X ~ g{x, mo, mi)) 



iy 



X — y + z [oj 

£pz = —[3y + F sin ur 

We choose the angular frequency tu = vr and define the strength by looking at 
the 1-d bifurcation diagram as a function of F. Corresponding to the windows 
in the diagram of Fig. 32, we consider following three types of drive control. 



Drive control-a : F = 1.4 
Drive control-b : F = 1.17 
Drive control-c : F = 1.7 



6.1 The drive control-a 




^.^ o.e i 
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Fig. 32. The 1-d bifurcation diagram as 
a function of the amplitude F of the 
drive control. 




Fig. 33. The drive control-a of the dou- 
ble scroll. 
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Fig. 34. The x coordinate of the trajectory of drive control-a type. 

A simulation data of F = 1.4, cj = vr is shown in Fig.33. The first derivative of 
the Poincare mapping is calculated from r = 143 to r = 543 where the initial 
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point is revisited. Its eigenvalues change from a pair of complex and a real to 
three real periodically. Similar behavior was observed [18] in the control of the 
Lorenz equation[19]. 
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Fig. 35. The real part of eigenvalues Fig. 36. The imaginary part of eigen- 
drive control-a type. values drive control-a type. 

In the interval from r = 338 to 382, the diffeomorphism is essentially in 2- 
dimensional plane and the two real eigenvalues satisfy the condition of the 

Newhouse region p < 1 < A < - and the orbit is absorped into the stable 

P 
periodic orbit. In the period of 387 < r < 405 products of the two eigenvalues 

pX are almost constant 0.91(1). 



6.2 The drive control-b 



A simulation data of F = 1.17, a; = vr is shown in Fig. 37. 
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Fig. 37. The drive control-b of the dou- 
ble scroll. 



Fig. 38. The x coordinate of the trajec- 
tory of drive control-b type. 



The first derivative of the Poincare mapping is calculated in the period from 
T = 1328 to r = 2728. Its eigenvalues change from a pair of complex and 
a real to three real periodically. The third eigenvalue is small and although 
the product of the complex eigenvalues become larger than 1, the absolute 



value of the product of three eigenvalues is small and the periodic trajectory 
is realized. Difference from the drive control-a is that the contraction mapping 
here is 3-dimensional. 



6.3 The drive control-c 
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Fig. 39. The drive control-c of the dou- 
ble scroll. 



Fig. 40. The x coordinate of the trajec- 
tory of drive control-c type. 



A simulation data oi F = 1.7,0; = tt is shown in Fig. 39. The first derivative 
of the Poincare mapping is calculated in the period between r = 1259 to 
r = 3259. One of its eigenvalues is consistent with zero. Interchanges of a pair 
of complex to real ones occur periodicaUy. The absolute value of the product 
of real eigenvalues pX exceeds 1 in the middle of the period but as it changes to 
a complex pair it is smaller than 1. Stretching in one direction and contracting 
in another direction occur when the eigenvalues are real and folding occurs 
when the eigenvalues are complex. Since the amplitude of the driving force is 
large, the oscillatory pattern of the complex eigenvalues is denser than that 
of the drive control-b type. 



7 Feedback control of the double scroll 



The application of feedback control is relatively simple. We choose the feed- 
back time delay to be Idr and simulate the following equation. 



dr 



X = a{y — X — g{x, nio, mi)) 



£y = x-y + z 

iz = -(3y + c(3y[T] 



(9) 



We select the strength c by looking at the window in the 1-d bifurcation 
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pattern as a function of c, as shown in Fig. 41. The trajectory obtained by 
using c = 0.1 is shown in Fig.42. This pattern is similar to the heterochnic 
orbit shown in [17]. 







Fig. 41. The 1-d bifurcation diagram 
as a function of the parameter c of the 
feedback control. 



Fig. 42. The feedback control of the 
double scroll. 
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Fig. 43. The x coordinate of the trajectory of feedback control type. 

The first derivative of the Poincare mapping is calculated from r = 64 to 
r = 2676 where the initial point is revisited. The x coordinate of the trajectory 
is shown in Fig. 43. The eigenvalues of the first derivative of the Poincare 
mapping changes from a pair of complex and a real to three real periodically. 
Their behavior is similar to that of the drive control-c type. 



8 Impulse control of the spiral 



The chaotic spiral can be obtained for the same parameters as those of double 
scroll written after eq.(2) except G = 0.65. The trajectory is shown in Fig. 44. 
We apply the impulse- 1 type whose strength is fixed by looking for windows 
in the 1-d bifurcation diagram as a function of the pulse strength shown in 
Fig. 45. Observing that a crisis occurs at Sy ~ —0.008, we apply the following 
control: 
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• Impulse-1 type: When x passes 1 from x„ > x„+i and y < —0.5 shift y to 
y - 0.0088. 

The Fig. 44 shows the chaotic pattern and the Fig.46 shows the result of 
impulse-1 type control. 

Although 2 pulses are triggered in the control of the double scroll, where hete- 
roclinic tangency occurs, only one pulse is triggered in the impulse-1 type con- 
trol, where homoclinic tangency occurs. The crossing of the orbit on W^{KJ) 
and that on iy*(A_) occurs at the fixed point P_ and the shift of the orbit on 
IV" (A_) into the region where countable periodic orbits exist makes the orbit 
ultimately tangential to W^{KJ) and the periodic orbit is realized. The con- 
trolled trajectory in the y — z plane and the homoclinic orbit started from the 
vicinity of the fixed point P'_ = (2.0 - 0.00001,0, -2.0) are shown in Fig.48. 
The X coordinate of the 3-cycle trajectory is shown in Fig.47. 




Fig. 44. The trajectory of the spiral 
without impulse control. 




Fig. 45. The 1-d bifurcation diagram as 
the function of the height of the im- 
pulse. 
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Fig. 46. The trajectory of spiral im- Fig. 47. The x coordinate of the trajec- 

pulse-1 type. tory of the impulse-1 type. 

The eigenvalues depend on the initial condition. But the presence of the 
stretching, contraction and rotation (bending) is the same. The mapping of 
the tetrahedron Af is almost parallel to the x — z plane and the shift of Az 
observed in the double scroll is absent. 
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Fig. 48. Tangency of the impulse-1 type controlled orbit and the homoclinic orbit. 



9 Discussion and outlook 



We showed that the Shil'nikov's chaos in the Chua's circuit can be controlled 
by 1) imposing impulse or 2) applying phase synchronization or 3) perform- 
ing feedback. The strength of the perturbation was defined by the position 
of windows in the 1-d bifurcation diagram. The perturbation on the unstable 
manifold in the window region does not destroy the structure of the stable 
manifold and the trajectory was absorbed into the sink. Due to homoclinicity, 
the orbit is expelled from the sink to an unstable manifold, but the pertur- 
bation again put the orbit to the sink and the periodicity is realized. In the 
control of spiral chaos, tangency of the homoclinic orbits predicted by the 
Newhouse's theorem was applied. In the control of the double scroll, hetero- 
clinic tangency accompanied by the shift in the z direction was utilized. 

This kind of control could be applied to the communication systems by defin- 
ing an appropriate coding function to the binary symbol sequence correspond- 
ing to rotations around P+ or -P-[l]. We observed that the number of cycles 
around the two fixed points depends on the detailed adjustment of the ca- 
pacitance and inductance. The shape of the pulse which is not a exact delta 
function modifies the destination of the orbit in the region near the fixed 
point. Various possible controlled double scroll due to heteroclinic tangency 
of Shil'nikov chaos and existence of several windows region in the 1-d bifurca- 
tion diagram allows encoding characters on each cyclic oscillation and apply 
to communication system. 

The chaos control method using the information of the 1-d bifurcation diagram[20] 
can be applied to other systems hke Rossler equation[22] and Langford equation. [21]. 
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